Shear-induced criticality near a liquid-solid transition of colloidal suspensions 



O 

(N 
-t— > 
O 

O 

S 1 

o 

B 

i 

-t— > 
o3 



o 
o 



d 
o 



X 



Masamichi J. Miyama 1 and Shin-ichi Sasa 1 

1 Department of Pure and Applied Sciences, University of Tokyo, Komaba, Meguro-ku, Tokyo 153-8902, Japan 

(Dated: October 8, 2010) 

We investigate colloidal suspensions under shear flow through numerical experiments. By measur- 
ing the time-correlation lunction of a bond-orientational order parameter, we find a divergent time 
scale near a transition point from a disordered fluid phase to an ordered fluid phase, where the order 
is characterized by a nonzero value of the bond-orientational order parameter. We also present a 
phase diagram in the (p, 7 CX ) plane, where p is the density of the colloidal particles and 7 ox is the 
shear rate of the solvent. The transition line in the phase diagram terminates at the equilibrium 
transition point, while a critical region near the transition line vanishes continuously as 7™ — > 0. 
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Introduction: A crystal phase is distinguished from a 
liquid phase by a translational and rotational symmetry 
breaking in space. Since there exists an order parameter 
associated with the symmetry-breaking, the nature of the 
transition to a crystal phase is rather different from that 
of gas-liquid transitions, although both are first-order 
transitions from the viewpoint of thermodynamics. In 
short, the transition to a crystal phase is classified as a 
symmetry-breaking first-order transition. 

The correlation time of fluctuations does not diverge 
near symmetry-breaking first-order transitions. This is 
in sharp contrast to the case of a symmetry-breaking 
second-order transition, which exhibits a divergent time 
scale, as is observed in the case of decreasing temperature 
around a critical point of liquid-gas transitions. On the 
basis of standard understanding, in the present Letter, 
we argue that a divergent time scale of steady-state fluc- 
tuations may appear near an equilibrium crystallization 
(symmetry-breaking first-order transition) point when a 
non-equilibrium condition is imposed on the system. 

Concretely, we study colloidal suspensions under shear 
flow through numerical experiments. Since the pioneer- 
ing work by Ackerson and Clark [l[ , there have been ex- 
tensive studies related to the crystallization of colloidal 
suspensions under shear flow [24a|- In particular, aphase 
diagram was obtained by numerical experiments [2| and 
laboratory experiments [3( , together with numerical real- 
izations of crystal-liquid coexistence under shear flow [1] . 
In order to clarify the microscopic mechanism for the 
transition to a crystal, the kinetics of homogeneous nu- 
cleation under shear flow was also investigated Q . How- 
ever, to our knowledge, shear-induced criticality near an 
equilibrium crystallization point has never been reported, 
except for our preliminary observation 

In the present study, we focus on time scales associ- 
ated with the relaxation from a bond-orientational or- 
dered state in a disordered regime. We first note that 
the relaxation time, even in equilibrium cases, diverges 
near the melting point. We characterize the parameter 
dependence of the relaxation time quantitatively by mea- 
surement of the time scries of a bond-orientational order 



parameter. The result is well fitted by the Vogel-Fulcher 
law, which suggests the existence of a nucleation process 
of disordered fluid regions in a crystal state. However, 
since crystals do not appear spontaneously in the dis- 
ordered phase, the divergent time scale has never been 
observed in steady state fluctuations. In contrast, in non- 
equilibrium systems under shear flow, the relaxation time 
exhibits a power-law divergence near a transition point 
to an ordered fluid, which suggests the existence of crit- 
ical slowing down. The power-law divergent time scale 
is also observed in steady state fluctuations. We refer to 
this novel phenomenon as shear-induced criticality near 
a liquid-solid transition. 

Model: We investigate N colloidal particles that are 
suspended in a solvent fluid confined to an L x L x L cu- 
bic box. We impose planar Couette flow on the solvent 
and choose the ir-axis and the z-axis to be the directions 
of the shear velocity and the velocity gradient, respec- 
tively. Concretely, the velocity profile of the solvent is 
assumed to be given as (7° x z,0,0). We impose periodic 
boundary conditions along the a:-axis and the y-axis and 
introduce two parallel walls so as to confine particles in 
the z direction. 

Let r,, i — 1, ••• ,N, be the position of particle i. 
The potential energy of particles U({fj}jLi) consists of 



two parts, Ei<j U i\ r i ~ r j\) and Eif 7 (n)- The 
former describes the interaction potential among parti- 
cles, where £/ LJ (r) = 4e((er/r) 12 - (cr/rf) - U cutoS for 
r < r c with cut-off length r c and J7 C utoff = 4e((cr/r c ) 12 — 
(cr/r c ) 6 ), while £/ LJ (r) = otherwise. [/ wall (r;) rep- 
resents the wall potential and is given by £/ wall (ri) = 
u WCA (r* - L/2 ± Zi) for L/2±z, < r* with the Weeks- 
Chandler- Andersen potential u WCA 0, while f7 wall (r,) = 
otherwise. 

We define momentum of the i-th particle relative to 
the shear flow as Pi(t) = mr^(t) — m A f ex Zi(t)e x , where 
to is the mass of a single particle. We then assume the 
equation of motion for the particles as 



dt 



dU({r 3 }f=i 

dn 



p 
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where & = (£f , , £?) represents thermal noise satisfying 

(£?(*)<f(*0) = 2CfcBT^<5 Q,3 (5(t-t'), is the Boltzmann 
constant, T is the temperature of the solvent, and C is the 
friction coefficient. The superscripts a and j3 represent 
Cartesian components. (See Ref. [§}.) 

In numerical simulations, all the quantities are con- 
verted to dimcnsionless forms by setting m = a = e = I . 
We fix T = 1.5, (/k B T = 1, N = 1024, r c = 2.5a, and 
r* = 0.5<7, and treat p and 7 ex as control parameters. We 
discretize ([I} according to the reversible system propaga- 
tor algorithm method [9( with time step At = 1/256. In 
the present Letter, (• • • ) represents the statistical average 
in steady states. 

Order parameters: We define an order parameter 



that characterizes a rotational symmetry breaking 10, 
ll| . Using the Delaunay triangular decomposition [12 1 
on a particle configuration, we determine neighbor- 
ing particles for a given particle i. The collection 
of edges that extend from in the Delaunay trian- 
gular is denoted by (rij)?=i i where ns(i) represents 
the number of neighbors of particle i. From this 
collection, we define a 13-dimensional vector q§(i) = 

(36,-6(0: • • • ) Q6,m{i), ■ ■ • 96,6 («)) &S 



96,m(0 



n B (i) 



(2) 



where Yq^ is the spherical harmonics function of degree 
six. Then, the bond-orientational order is qualified by 



q«,, 



1 N 



(i). 



(3) 



Here, (qe, m ) = if the rotational symmetry is not broken, 
while (q~6,m) 7^ in the thermodynamic limit ./V — > 00 
when the bond-orientational order emerges. In order to 
detect symmetry-breaking, it is convenient to measure 
the magnitude of the vector QQ^rn- 

Following the standard 

convention, we define 




1/2, 



(?6,m(?6,m 



(4) 



Note that Qq ~ 0(X/s/N) in the disordered phase, while 
Qq ~ O(l) in the ordered phase, when N —> 00. 

On the left-hand side of Fig. Q] we show Qq as a func- 
tion of /9 for several values of 7 CX . The figure indicates the 
existence of an ordered state with Qq ~ 0(1) in a high- 
density regime for each 7 CX . In particular, the transition 
to the ordered phase is quite sharp when 7™ = 0, while 
the transition width becomes wider as -y cx is increased. 
For a tentative value of the transition point, we define 
p q as the density such that Qq = 0.2. We display p q as 
a function of 7 ex in the right-hand side of Fig. [TJ More 
precise determination of the functional forms of Qq will 
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FIG. 1: (Color online) Left: Q§ as a function of p for several 
values of 7 CX . Right: p q as a function of 7 CX . 



be obtained by investigating larger systems. Note that 
in the thermodynamic limit without shear flow, non-zero 
Qe emerges continuously for a density at which a crys- 
tal can coexist with a liquid. When we ignore the co- 
existence phase, Qq exhibits a discontinuous transition, 
which is observed for the system under constant pres- 
sure. In the present Letter, putting aside phenomena 
associated with the coexistence phase, we focus on the 
question of how the nature of the symmetry-breaking 
discontinuous-transition is modified by the influence of 
shear flow. 

In the equilibrium case, the ordered phase corresponds 
to a crystal. However, in the non-equilibrium cases, since 
the shear flow drives particles, particles may flow even in 
the ordered phase with Qq ^ 0. We then measure the 
x-component of the velocity averaged over a region with 
an interval [z + 0.5, z — 0.5] in the z direction, which is 
denoted by v(z). Examples of v{z) for several values of 
7° x with p = 1.1 fixed are shown in the inset of Fig. [2] 
We then determine the shear rate 7 of particles by fitting 
the slope of the velocity profile v(z) in the region —1.5 < 
z < 1.5. The obtained shear rates 7 are plotted for j CK 
in Fig. [5] Although the flow might cease at some value 
of 7° x , the determination as to whether the cross-over is 



actually singular is a delicate problem. (See Ref. [14 1 
for a related discussion.) For any case, the cross-over 
points are located at a higher density than p q when 7° x > 
0. (See Fig. O) Therefore, as we are concerned with 
behaviors near the transition point at which the order 
parameter Qq ~ O(l) appears, we may assume that an 
ordered fluid is observed in the ordered phase in the non- 
equilibrium cases. 

Transition to the ordered fluid: Next, we character- 
ize the nature of the transition to the ordered phase. 
First, in order to observe a clear difference between the 
equilibrium and non-equilibrium cases, we measure the 
relaxation time r rc i at which Qq reaches a value of 0.05, 
starting from a crystal state, which Qq s» 0.35. Note 
that r re i can be measured only in the disordered phase. 
In Fig. [H we show t ic \ as a function of p for several val- 
ues of 7 ex , where we set the maximum waiting time to 
t = 10000. The results indicate the existence of a charac- 
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FIG. 2: (Color online) Shear rate of particle flow as a function 
of 7 CX for several values of p. Each point is obtained from 
the velocity profiles v(z). Examples of velocity profiles for 
different 7° x with p — 1.1 fixed are shown in the inset. 
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FIG. 3: (Color online) Phase diagram in the (p,7° x ) plane. 
The green and cross symbols represent the transition line p = 
Pq(7 ex ) between the disordered fluid phase and the ordered 
fluid phase. The parameter values for realizing 7 = 10~ 3 , 
below which flow appears to cease, are also plotted as red 
plus symbols. 



teristic density pd at which the relaxation time diverges 
for each value of 7 ex . 

Let us determine the functional form of T ro i with the 
value of pd- First, as shown in the inset of Fig. [5j r re i for 
the equilibrium case is well fitted by the Vogel-Fulcher 
law 



Trel — T exp 



.4 



Pd 



(5) 



A phenomenological argument may be developed for the 
nucleation of a disordered domain, by which ([5]) may be 
understood. (See for example Ref. 13J for a demonstra- 
tion of a g-states Potts model.) In contrast, Fig. [5] indi- 
cates that T ro i for the systems under shear flow follows a 
power-law form 



Trol 



5(7 ex )(Pd(7 cx )-pr 



(0) 



where £ ~ 1.6. This suggests that the divergent behavior 
does not originate from the nucleation of disordered re- 
gions but may be related to critical slowing down. Note 
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FIG. 4: (Color online) Relaxation time r re i as a function of 
density p for several values of 7 ex . 
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FIG. 5: (Color online) Fitting of functional forms of r rc i. r ro i 
versus pa — p with a log-log plot for -y cx = 0.1, 0.316, and 
1.0. pd is a fitting parameter, the value of which is estimated 
as 1.006, 1.037, and 1.121, respectively. The guide line rep- 
resents ([6]) with = 1.6. Inset: logr rc i versus pd — p with a 
log-log plot for 7° x = 0. The guide line represents © with 
A = 0.1 and p d = 1.003. 



that the pre-factor B in © depends slightly on 7 CX in 
the form B ~ (^ycx^o.3^ ag s j 10wn m ^ e j nse t f Fig. [B] 
Let p w (7 ex ) be a typical width of the power-law region 
for 7° x . Then, on the basis of the dimensional analysis, 
it is expected that 



7 ex 7Vei 



Pd(7 cx ) ~P 



< 



(7) 



By assuming p w ~ (7 CX ) X in ([7]), we obtain B ~ 
^ eX )xC-i ) w hich leads to x - 1-3/1.6 > 0. This means 
that a critical region for the system with finite 7 CX be- 
comes narrower for smaller 7° x and vanishes in the equi- 
librium system. See the schematic phase diagram in 
Fig. [6] 

We now note that such a divergent time scale is never 
observed in the stationary state of the equilibrium sys- 
tem. In order to confirm this explicitly, we measure the 
time correlation function defined by 



C(t) = 



Ei 

m— — 6 



(g6,m(M<Z6,m(*0 + t)) - (q 6 ,; 



(8) 
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FIG. 6: (Color online) Schematic phase diagram with the 
critical region in the (p-7 cx ) plane. Inset: _B(7 CX ) as a function 
of 7 CX . The guide line represents a power- law function with 
exponent 0.3. 
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FIG. 7: (Color online) Correlation time r corr in steady states 
as a function of p. The melting time r rc i is also superimposed 
for comparison. (Left:) -y cx = 0, and (Right:) -y cx = 1.0. 



We determine the correlation time r c from the fitting of 
the exponential decay rate of C(t). On the left-hand side 
of Fig. r c is shown as a function of p in the disordered 
regime. Indeed, the correlation time does not diverge. 
For reference, we superimpose the data of the relaxation 
time 7" re i. In contrast to the equilibrium case, as shown in 
the right-hand side of Fig. the correlation time in the 
system under shear flow diverges in a manner similar to 
T Ie i. These results indicate that the symmetry-breaking 
transition to the ordered fluid accompanies a critical phe- 
nomenon. This is the main claim of the present Letter. 

Concluding remarks: Before ending this Letter, we 
address three considerations. First, we conjecture that 
the fluctuation of q 6 possesses a divergent length scale. 
By investigating the manner of divergences of several 
quantities for systems of different sizes, the universality 
class for this phenomenon may be determined. 

Second, with regard to the universality problem, we 
are also interested in a simple mathematical model in the 
same universality class. For example, it might be possible 
to propose a model describing a stochastic time evolution 



of the coarse-grained order parameter field. The first 
problem in this direction is to derive the value of £ using 
a phenomenological argument. 

Finally, in all of the arguments presented above, the 
coexistence phase is ignored. In order to extract more 
precise results, it might be better to investigate systems 
under constant pressure. A study of such systems of 
larger sizes will be performed in the future. 

In summary, we have investigated colloidal suspensions 
under shear flow. We have found that the critical tran- 
sition line starts from the liquid-solid transition point in 
the equilibrium system without shear flow. This novel 
phenomenon, referred to as shear-induced criticality, will 
be investigated from several viewpoints. 
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